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THREE-DIMENSIONAL SPATIAL FILTERING APPARATUS AND METHOD 

BACKGROUND OF THE INVENTION 

The present invention relates to a three-dimensional spatial filtering 
5 apparatus and method, and more particularly to a three-dimensional spatial 
filtering apparatus and method for conducting three-dimensional spatial filtering 
on a three-dimensional image. 

Three-dimensional spatial filtering is conducted on a three-dimensional 
image for improvement of image quality of the three-dimensional image, for 
10 example. Several types of the three-dimensional spatial filtering include, for 
example, low-pass filtering, edge enhancement, and correlation detecting 
filtering. The three-dimensional spatial filtering is achieved by scanning image 
data by a three-dimensional spatial filter and superposing them on each other. 
The superposition is sometimes referred to as convolution. Weighting factors 
15 for the convolution are sometimes referred to as filter coefficients. 

Conventionally, the filter coefficients for the three-dimensional spatial 
filtering are fixed at predetermined values. Thus, low-pass filtering, for 
example, is applied at a uniform strength throughout a three-dimensional image. 
When the three-dimensional image is comprised of a plurality of 
20 tomographic images obtained by, for example, an X-ray CT (computed 
tomography) apparatus, then in the case where the tomographic images are of 
low-dose X-rays, noise contained in the images may be larger in a portion having 
a larger pixel value. This results from the fact that the pixel value is represented 
by a CT number. Specifically, since the CT number is a numerical value 
25 indicating attenuation of X-rays, a larger pixel value represents weaker X-rays 
having a poorer S/N (signal-to-noise ratio) and increased noise. 

When, for example, conducting low-pass filtering as three-dimensional 
spatial filtering with fixed filter coefficients on a three-dimensional image of the 
foregoing nature, filtering suitable for a portion having a larger pixel value and a 
30 poorer S/N is excessive for a portion having a smaller pixel value and a better 
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S/N. On the other hand, filtering for a portion having a smaller pixel value and 
a better S/N is insufficient for a portion having a larger pixel value and a poorer 
S/N. Such a case is one in which suitable filtering cannot be achieved 
throughout a three-dimensional image. 

5 

SUMMARY OF THE INVENTION 

It is therefore an object of the present invention to provide a 
three-dimensional spatial filtering apparatus and method for suitably conducting 
three-dimensional spatial filtering throughout a three-dimensional image even if 
10 it has a varying property dependent upon a pixel value. 

(1) The present invention, in accordance with one aspect for solving the 
aforementioned problem, is a three-dimensional spatial filtering apparatus for 
conducting filtering on a three-dimensional image by a three-dimensional spatial 
filter, characterized in comprising: coefficient adapting means for changing 

15 coefficients of said three-dimensional spatial filter according to a property of a 
pixel value in said three-dimensional image. 

(2) The present invention, in accordance with another aspect for solving 
the aforementioned problem, is a three-dimensional spatial filtering method for 
conducting filtering on a three-dimensional image by a three-dimensional spatial 

20 filter, characterized in comprising: adapting coefficients of said 

three-dimensional spatial filter according to a property of a pixel value in said 

three-dimensional image. 

In the invention of the aspects described in (1) and (2), since coefficients of 

a three-dimensional spatial filter are changed according to a property of a pixel 
25 value in a three-dimensional image, the three-dimensional spatial filtering can be 

suitably achieved throughout a three-dimensional image. 

Preferably, said property is a pixel value of a pixel-of-interest in the 

three-dimensional image convoluted with the filter, insofar as the value of the 

pixel-of-interest is to take effect on the filter coefficients. 
30 Preferably, said property is an average value of pixel values of the 
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pixel-of-interest in the three-dimensional image convoluted with the filter and its 
neighboring pixels, insofar as the average value of pixel values of the 
pixel-of-interest and its neighboring pixels is to take effect on the filter 
coefficients. 

5 Preferably, said property is a maximum value of pixel values of the 

pixel-of-interest in the three-dimensional image convoluted with the filter and its 
neighboring pixels, insofar as the maximum value of pixel values of the 
pixel-of-interest and its neighboring pixels is to take effect on the filter 
coefficients. 

10 Preferably, said property is a rninimum value of pixel values of the 

pixel-of-interest in the triree-dimensional image convoluted with the filter and its 
neighboring pixels, insofar as the rninimum value of pixel values of the 
pixel-of-interest and its neighboring pixels is to take effect on the filter 
coefficients. 

15 Preferably, said property is a median value of pixel values of the 

pixel-of-interest in the three-dimensional image convoluted with the filter and its 
neighboring pixels, insofar as the median value of pixel values of the 
pixel-of-interest and its neighboring pixels is to take effect on the filter 
coefficients. 

20 Preferably, said property is a maximum value of absolute values of 

differences between a pixel value of the pixel-of-interest in the three-dimensional 
image convoluted with the filter and pixel values of its neighboring pixels, 
insofar as the maximum value of absolute values of differences between pixel 
values of the pixel-of-interest and its neighboring pixels is to take effect on the 

25 filter coefficients. 

Preferably, said property is a maximum value of squares of differences 
between a pixel value of the pixel-of-interest in the three-dimensional image 
convoluted with the filter and pixel values of its neighboring pixels, insofar as the 
maximum value of squares of differences between pixel values of the 

30 pixel-of-interest and its neighboring pixels is to take effect on the filter 
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coefficients. 

Preferably, said property is a minimum value of absolute values of 
differences between a pixel value of the pixel-of-interest in the three-dimensional 
image convoluted with the filter and pixel values of its neighboring pixels, 
5 insofar as the minimum value of absolute values of differences between pixel 
values of the pixel-of-interest and its neighboring pixels is to take effect on the 
filter coefficients. 

Preferably, said property is a minimum value of squares of differences 
between a pixel value of the pixel-of-interest in the three-dimensional image 
1 0 convoluted with the filter and pixel values of its neighboring pixels, insofar as the 
minimum value of squares of differences between pixel values of the 
pixel-of-interest and its neighboring pixels is to take effect on the filter 
coefficients. 

Preferably, said property is a median value of absolute values of 
15 differences between a pixel value of the pixel-of-interest in the three-dimensional 
image convoluted with the filter and pixel values of its neighboring pixels, 
insofar as the median value of absolute values of differences between pixel 
values of the pixel-of-interest and its neighboring pixels is to take effect on the 
filter coefficients. 

20 Preferably, said property is a median value of squares of differences 

between a pixel value of the pixel-of-interest in the three-dimensional image 
convoluted with the filter and pixel values of its neighboring pixels, insofar as the 
median value of squares of differences between pixel values of the 
pixel-of-interest and its neighboring pixels is to take effect on the filter 

25 coefficients. 

Preferably, the coefficients of the three-dimensional spatial filter are 
adapted depending upon, as said property, a property of a pixel value of the 
pixel-of-interest in the three-dimensional image convoluted with the filter, 
insofar as the property of the pixel value is to take effect on the filter coefficients. 
30 It may sometimes be preferable from the point of reflecting the property of 
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the pixel value in the filter coefficients to employ the standard deviation of the 
pixel-of-interest in the three-dimensional image convoluted with the filter and its 
neighboring pixels as one indicator of said property. 

It may sometimes be preferable from the point of reflecting the properties 
5 of the pixel-of-interest and its neighboring pixels in the filter coefficients to 
employ a sum of absolute values of differences between a pixel value of the 
pixel-of-interest in the three-dimensional image convoluted with the filter and 
pixel values of its neighboring pixels as an indicator of said property. 

It may sometimes be preferable from the point of reflecting the properties 
10 of the pixel-of-interest and its neighboring pixels in the filter coefficients to 
employ a square sum of differences between a pixel value of the pixel-of-interest 
in the three-dimensional image convoluted with the filter and pixel values of its 
neighboring pixels as an indicator of said property. 

Therefore, the present invention can provide a three-dimensional spatial 
15 filtering apparatus and method for conducting three-dimensional spatial filtering 
with coefficients dependent upon a pixel value. 

Further objects and advantages of the present invention will be apparent 
from the following description of the preferred embodiments of the invention as 
illustrated in the accompanying drawings. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a block diagram of an exemplary apparatus in accordance with 
one embodiment of the present invention. 

Figure 2 is a block diagram of an image processing section. 
25 Figure 3 is a conceptual diagram of three-dimensional image data. 

Figure 4 is a conceptual diagram of a three-dimensional spatial filter. 
Figure 5 is a conceptual diagram of a convolution filter. 
Figure 6 is a diagram showing an exemplary coefficient setting in the 
three-dimensional spatial filter. 
30 Figure 7 is a diagram showing an exemplary coefficient setting in the 
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three-dimensional spatial filter. 

Figure 8 is a flow chart of an operation of the exemplary apparatus in 
accordance with one embodiment of the present invention. 

Figure 9 is a diagram showing an exemplary coefficient setting in the 
5 three-dimensional spatial filter. 

Figure 10 is a diagram showing an exemplary coefficient setting in the 
three-dimensional spatial filter. 

Figure 1.1 is a diagram showing an exemplary coefficient setting in die 
three-dimensional spatial filter. 
10 Figure 12 is a diagram showing an exemplary coefficient setting in the 

three-dimensional spatial filter. 

DETAILED DESCRIPTION OF THE INVENTION 

Embodiments of the present invention will now be described in detail with 

15 reference to the accompanying drawings. Figure 1 shows a block diagram of a 
three-dimensional spatial filtering apparatus. The apparatus is one 
embodiment of the present invention. The configuration of the apparatus 
represents an embodiment of the apparatus in accordance with the present 
invention. The operation of the apparatus represents an embodiment of the 

20 method in accordance with the present invention. 

As shown in Figure 1, the present apparatus has an image acquiring section 
102. The image acquiring section 102 is for acquiring a three-dimensional image 
of a certain subject. The image acquiring section 102 employed is, for example, 
an X-ray CT apparatus. The X-ray CT apparatus captures a plurality of 

25 tomographic images of the subject by a helical scan or the like, that is, captures a 
three-dimensional image. 

Instead of an X-ray CT apparatus, the image acquiring section 102 may be a 
magnetic resonance imaging (MRI) apparatus or an ultrasonic imaging apparatus, 
or a three-dimensional imaging apparatus for a non-medical application such as 

30 non-destructive inspection or physical phenomenon analysis, e.g., an industrial 
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X-ray CT apparatus or laminography apparatus. 

The three-dimensional image acquired by the image acquiring section 102 
is input to an image processing section 104. The image processing section 104 
conducts three-dimensional spatial filtering on the input three-dimensional 
5 image. The image processing section 104 is connected with a display section 106 
and an operating section 108 so that a user can interactively operate the image 
processing section 104 therethrough as needed. 

The image processing section 104 is comprised of, for example, a computer. 
Figure 2 shows an exemplary configuration of the image processing section 104 
10 as a block diagram. As shown, the image processing section 104 is comprised of 
a CPU (central processing unit) 402, a main memory 404, an external memory 406, 
an image input interface 408, a display output interface 410, and an operation 
input interface 412, all connected with one another via a bus 414. 

The image input interface 408 is connected with the image acquiring 
15 section 102. The display output interface 410 is connected with the display 
section 106. The operation input interface 412 is connected with the operating 
section 108. 

The external memory 406 stores a program for operating the CPU 402. 
The CPU 402 reads the program into the main memory 404 and executes it to 
20 achieve predefined functions as the image processing section 104. 

The external memory 406 also stores three-dimensional image data 
gathered from the image acquiring section 102 via the image input interface 408. 
Figure 3 is a conceptual diagram of the three-dimensional image data stored in 
the external memory 406. As shown, the three-dimensional image data forms a 
25 three-dimensional pixel matrix. 

When directions of three axes of the three-dimensional pixel matrix are 
represented as x, y and z, the matrix size is, for example, 512 pixels both in the x- 
and y-directions and 100 pixels in the z-direction. That is, it is a 512 x 512 x 100 
matrix. The matrix corresponds to a set of 100 consecutive slices of tomographic 
30 images each having 512 x 512 pixels. 
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The three-dimensional image is subjected to three-dimensional spatial 
filtering by the CPU 402. The three-dimensional spatial filtering is carried out 
by convoluting the image data with a three-dimensional spatial filter. 

Figure 4 shows an example of the three-dimensional spatial filter. As 
5 shown, the three-dimensional spatial filter has coefficients arranged in a 3 x 3 x 3 
configuration. It should be noted that the size of the three-dimensional spatial 
filter is not limited to 3 x 3 x 3 but may be an appropriate one such as 5 x 5 x 5 or 
7x7x7, or 2x2x2 or 4x4x4. Moreover, instead of the size being equilateral 
in the x-, y- and z-directions, it' can be unequilateral. 
10 Filter coefficients are defined for the three-dimensional spatial filter. The 

three-dimensional spatial filter is sometimes referred to as a convolution filter. 
Figure 5 shows an example of the convolution filter. 

As shown, the convolution filter performs a convolution calculation in a 
one-to-one correspondence with each pixel in a local region having the same size. 
15 At each pixel position, a coefficient is disposed. Such convolution filter 
coefficients are stored in the external memory 406 beforehand. 

The three-dimensional spatial filtering is an operation of scanning an image 
with the convolution filter and calculating a value by a product-sum operation 
on the convolution filter and the pixel values of the corresponding pixels. Thus, 
20 an image of results of the convolution product-sum operation by the convolution 
filter is obtained. This operation is the three-dimensional spatial filtering. 

The characteristic of the filtering is determined by the coefficient values in 
the three-dimensional spatial filter and their positions. Figure 6 shows an 
example of the coefficient values in the three-dimensional spatial filtering and 
25 their positions. 

Figure 6(a) shows a case in which identical coefficients, 1/27, are disposed 
in all the pixel positions. Such a coefficient setting achieves averaging, i.e., 
low-pass filtering. Figure 6(b) shows a case in which coefficients at pixel 
positions along the y-direction in the center are all defined as zero, coefficients at 
30 pixel positions along the y-direction on one side are all defined as minus one, and 
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those on the other side are all defined as one. This achieves first-order 
differentiation, i.e., edge enhancement, in the x-direction. Figure 6(c) shows a 
case in which a central pixel is defined as zero, coefficients neighboring it in the 
x-, y- and z-directions are defined as minus one, and other pixels are defined as 
5 zero. This achieves second-order differentiation. 

In the three-dimensional spatial filtering, it is also possible to set 
coefficients for effecting edge enhancement in a direction oblique with respect to 
all of the x- y- and z-axes. Figure 7 shows an example of such a coefficient 
setting. Figures 7(a), (b) and (c) show an exemplary coefficient setting of three 
10 layers of x-y spatial filter coefficients, respectively, arranged in the z-axis in the 
three-dimensional spatial filtering. 

In the first layer (z - 1), as shown in (a), coefficients of one are positioned 
along a diagonal direction declining toward the right, coefficients of one are 
positioned in an entire right portion with respect to the diagonal direction, 
15 coefficients of zero are positioned in a left portion except a corner, and a 
coefficient of minus one is positioned at the corner. 

In the second layer (z), as shown in (b), coefficients of zero are positioned 
along a diagonal direction declining toward the right, coefficients of one are 
positioned in an entire right portion with respect to the diagonal direction, and 
20 coefficients of minus one are positioned in an entire left portion. 

In the third layer (z + 1), as shown in (c), coefficients of minus one are 
positioned along a diagonal direction declining toward the right, coefficients of 
minus one are positioned in an entire left portion with respect to the diagonal 
direction, coefficients of zero are positioned in a right portion except a corner, 
25 and a coefficient one is positioned at the corner. 

By such coefficient arrangement, an edge lying in a diagonal direction 
declining toward the "lower-right" and in a direction obliquely crossing the 
three-dimensional spatial filter from the first layer to the third layer can be 
enhanced. 

30 The image processing section 104 conducts filtering while changing the 
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coefficient setting for the three-dimensional spatial filter according to the 
property of the pixel value or the pixel value itself of a pixel-of-interest. The 
property of the pixel value for the three-dimensional spatial filter includes the 
pixel value of the pixel-of-interest for the three-dimensional spatial filter. Other 
5 properties include the average value, maximum value, mirumum value, median 
value, and the like, of the pixel-of-interest and its neighboring pixels. 

Moreover, the property of the pixel value may be a value given by each of 
the following equations: 



o 



max\G(x, y, z) - G(x -i,y- j, z - k% 



10 




(1) 



max {G(x 9 y, z) - G{x - 1 



i,y-j\z- k)} 2 ; 



where 



*=i,-l; 

j = 1, - 1; and 



(2) 



min|G(x ? y, z) - G(x -i,y- j\ z - k)\ 




(3) 



15 
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min {G(x, y, z)-G(x- i, y-j,z- k)} 2 ; 
where 

/=i,-i; (4) 

7=1,- 1; and 



mod|G(x, y, z)-G(x- i, y-j,z- k% 
where 

» =1,-1; (5) 

y =1,-1; 

£ = 1, - 1, and 



mod {G(x 5 j/, z) - G(x - /, j/ -y, z - A)} 2 ; 

i" (6) 

j = 1,-1; and 
k = h-l. 



Equation (1) is an embodiment representing the maximum value of the 
absolute values of the differences between pixel values of the pixel-of-interest for 
the three-dimensional spatial filtering and other pixels. 
10 Equation (2) represents the maximum value of the squares of the 

differences between pixel values of the pixel-of-interest for the three-dimensional 
spatial filtering and other pixels. 

Equation (3) represents the minimum value of the absolute values of the 
differences between pixel values of the pixel-of-interest for the three-dimensional 
1 5 spatial filtering and other pixels. 

Equation (4) represents the minimum value of the squares of the 
differences between pixel values of the pixel-of-interest for the three-dimensional 
spatial filtering and other pixels. 
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Equation (5) represents the median value of the absolute values of the 
differences between pixel values of the pixel-of-interest for the three-dimensional 
spatial filtering and other pixels. 

Equation (6) represents the median value of the squares of the differences 
5 between pixel values of the pixel-of-interest for the three-dimensional spatial 
filtering and other pixels. 

Moreover, another example of the property of the pixel value in the 
three-dimensional spatial filtering may be represented by the standard deviation. 
Exemplary indicators of the property other than the standard deviation include a 
1 0 value given by each of the following equations: 



£ \G(x, y,z)-G(x- U y-j 9 z- k% 
where 

/ = (7) 
j = 1,-1; and 
k = 1, - 1, and 



X y* z ) - G ( x - '> y -J> z - k )} 2 > 

where 

i=i,-i; (8) 

j = 1, - 1; and 
A: = l,-1. 

15 

Equation (7) represents the sum of the absolute values of the differences 
between pixel values of the pixel-of-interest for the three-dimensional spatial 
filtering and other pixels. 

Equation (8) represents the square sum of the differences between pixel 
20 values of the pixel-of-interest for the three-dimensional spatial filtering and other 
pixels. 

Figure 8 shows a flow chart of an operation of the image processing section 
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104. The operation of the image processing section 104 is substantially an 
operation of the CPU 402. As shown, the three-dimensional spatial filter starts 
scanning. At Step 702, initialization is first conducted. By this step, the initial 
position of the pixel-of-interest is defined as x = 0, y = 0 and z = 0. 
5 Next, at Step 704, pixel values of the pixel-of-interest and its neighboring 

pixels are read from the three-dimensional image. 

Next, at Step 706, property determination is conducted. The property to 
be employed is determined beforehand as any one of the pixel value of the 
pixel-of-interest, maximum value of neighborhood, minimum value of 

10 neighborhood and the like, as described above. 

For example, if the property specified by the user is the pixel value of the 
pixel-of-interest for the three-dimensional spatial filter, a value G0(x, y, z) in the 
input image is read. Although the following description will be made on a case 
in which the pixel value of the pixel-of-interest for the three-dimensional spatial 

15 filter is employed as the property, the same applies to other cases. 

Next, at Step 708, coefficient setting for the three-dimensional spatial filter 
is conducted. The CPU 402 is an embodiment of the coefficient adapting means 
of the present invention. The coefficient setting for the three-dimensional 
spatial filter is achieved according to the property determined at Step 706, i.e., 

20 according to the pixel value of the pixel-of-interest, for example. The pixel 
value of the pixel-of-interest will be sometimes referred to as the pixel-of-interest 
value hereinbelow. 

The coefficient setting is achieved in the following manner, for example. 
When the three-dimensional image is comprised of a plurality of tomographic 

25 images acquired by an X-ray CT apparatus and the three-dimensional spatial 
filter is a low-pass filter, the coefficients are set so that the strength of filtering is 
greater for the larger pixel-of-interest value and is smaller for the smaller 
pixel-of-interest value. In other words, the coefficients are set so that averaging 
is strengthened for the larger pixel-of-interest value and is weakened for the 

30 smaller pixel-of-interest value. 
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Figure 9 shows an exemplary coefficient setting for the strong low-pass 
filtering. As shown, all the pixel positions in the three-dimensional spatial filter 
are assigned a coefficient of 1/27. By such a coefficient setting having the sum 
of the coefficients normalized to "1", pixel values of the pixel-of-interest and its 
5 neighboring pixels all convoluted by the same weight are obtained. This simple 
averaging corresponds to the stronger low-pass filtering. 

Figure 10 shows an exemplary coefficient setting for the weak low-pass 
filtering. As shown, the central pixel position and vertically and horizontally 
neighboring pixel positions in each layer of the x-y plane of the 

10 three-dimensional spatial filter are assigned a coefficient of 1/15, and pixels in 
the remaining portion are assigned a coefficient of zero. By such a coefficient 
setting having the sum of the coefficients normalized to "1", the degree of 
averaging is reduced by an amount corresponding to pixel values in a portion 
having the coefficient of zero that are excluded from the convolution. 

15 The larger the number of pixel positions assigned a coefficient of zero, the 

smaller the degree of averaging; and as shown in Figures 11(a) and (c), when 
only the central pixel is assigned a coefficient of 1/7 and the other pixels are 
assigned zero in the first and third layers, yet weaker averaging can be achieved. 
Moreover, when only the pixel-of-interest is assigned a coefficient of one and the 

20 remaining pixels are all assigned zero as shown in Figure 12, no change is 
effected. 

The setting of coefficients according to the pixel-of-interest value is 
conducted by a predetermined algorithm or calculation of an equation. 
Alternatively, it is possible to store beforehand a table representing 
25 correspondence of the pixel-of-interest value and the coefficient and conduct the 
coefficient setting with reference to the table. 

Next, at Step 710, convolution is conducted. The convolution is 
processing of conducting convolution on a pixel scanned by the 
three-dimensional spatial filter having the coefficients defined as described above, 
30 with such a three-dimensional spatial filter. The result of the convolution of the 
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pixel-of-interest and its neighboring pixels with the three-dimensional spatial 
filter gives a three-dimensional spatial filtering output value for the 
pixel-of-interest. 

Tomographic images captured by an X-ray CT apparatus have larger noise 
5 in a portion of a larger pixel value and smaller noise in a portion of a smaller 
pixel value as described earlier; and therefore, by strengthening smoothing for a 
larger pixel-of-interest value and weakening smoothing for a smaller 
pixel-of-interest value as described above, low-pass filtering for a plurality of 
tomographic images can be suitably achieved throughout the images. 

10 Unlike the tomographic images captured by the X-ray CT apparatus, a 

three-dimensional image such as that captured by an MRI apparatus, for example, 
which has a smaller noise in a portion of a larger pixel value and larger noise in a 
portion of a smaller pixel value, is subjected to, contrary to the above, weaker 
smoothing in a portion of a larger pixel-of-interest value and stronger smoothing 

15 in a portion of a smaller pixel-of-interest value, whereby low-pass filtering can be 
suitably achieved throughout the image. 

In first- or second-order differentiation, for example, in addition to the 
low-pass filtering, suitable processing can be achieved throughout a 
three-dimensional image by setting coefficients as described above. 

20 Next, at Step 712, a pixel value is written into an output image. The pixel 

value written constitutes a three-dimensional filtered output image. 

Next, at Step 714, x is set as x = x + 1. This moves the position of the 
pixel-of-interest by one in the x-direction. Next, at Step 716, a decision is made 
as to whether x = N. N is the number of pixels in the input image in the 

25 x-direction. If the decision results in No, the process goes back to Step 704 and 
sequentially conducts processing of Steps 704 — 716. This process is repeated 
until x = N is reached. 

When x = N is reached, x and y are set as x = 0 and y = y + 1 at Step 718. 
This moves the position of the pixel-of-interest by one in the y-direction. Next, 

30 at Step 720, a decision is made as to whether y = N. N is the number of pixels in 
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the input image in the y-direction. 

If the decision results in No, the process goes back to Step 704 and 
processing at Steps 704 — 716 are repeated until x = N is reached, y is 
incremented by one each time x = N is reached, and the process is repeated until 
5 y = N is reached. 

When y = N is reached, x, y and z are set as x = 0, y = 0, and z = z + 1 at 
Step 722... This moves the position of the pixel-of-interest by one in the 
z-direction. Next, at Step 724, a decision is made as to whether z = M. M is the 
number of pixels in the input image in the z-direction. 

10 If the decision results in No, the process goes back to Step 704 and 

processing at Steps 704 — 716 are repeated until x = N is reached, y is 
incremented by one each time x = N is reached, the process is repeated until y = 
N is reached, z is incremented by one each time y = N is reached, and the process 
is repeated until z = M is reached. When z = M is reached, the 

1 5 three-dimensional spatial filtering is completed. 

Many widely different embodiments of the invention may be configured 
without departing from the spirit and the scope of the present invention. It 
should be understood that the present invention is not limited to the specific 
embodiments described in the specification, except as defined in the appended 

20 claims. 



